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ON THE EXTENSION OF A THEOREM OF POINCARfi FOR 
DIFFERENCE-EQUATIONS * 

BY 

EDWARD B. VAN VLECK 

One of the most useful theorems relating to difference equations is the well- 
known theorem of PoiNCAKEf for the equation 

(1) u(n+r) + Pi{n)u{n+r- 1)+ \- Pr (n) u (n) = , 

in which the coefficients Pi (ra) possess limits Ai for n = -\- co. Poincake 
showed, namely, that if the roots of the equation 

(2) y^ + A^y^'+---+Ar=0 

are no two of equal modulus, the ratio u (n) / u {n — 1) has a limit for 
w = + 00 , and this limit is a root of (2), in general the root of greatest 
modulus.^ 

* Presented to the Society at the Chicago meeting, April 29, 1911. 

tAmerican Journal of Mathematics, vol. 7 (1885), p. 213. For careful anal- 
ysis and further development of the theorem, see Perron, Crelle's Journal, vol. 136 
(1909), p. 17 and vol. 137, p. 6. 

X In 1904, in an article on continued fractions in these Transactions, p. 253, I gave a 
proof of Poincare's theorem for r = 2 which was suitable for application to the class of con- 
tinued fractions there under consideration. This proof, and consequently its immediate appli- 
cability, was not fully understood by Phingsheim and called forth an interesting memoir from 
him (Sitzungsberichte der Bayerischen Akademie, math.-physika- 
lischer Klasse, 1910, p. 3) devoted in large part to rederiving my results. Later 
after oral and written communication between us he has recognized the legitimacy of my 
proof (Ibid., 1911, p. 64). My argument was somewhat condensed, but not more so, it seemed 
to me, than is current in articles of an investigational character. A few additional lines, which 
it is worth while to add now, would have removed any difficulty. 

When Pringsheim reaches the following theorem (p. 46), he inadvertently fails to at- 
tribute it to me, though the attribution is intended in the introduction : 

// in the continued fraction 

^(^>- 1+ 1 + 1 + 

we have Km. a„ = fc for n = a:i , the continued fraction will converge over the entire plane except 
(1) along the whole or a part of a rectilinear cut drawn from z = — i k to z = oo with an argument 
equal to that of the vector from the origin to z = — \k, and except possibly (2) at certain isolated 
points pi, p2, • • • . Within the plane so cut the continued fraction is holomorphic except at the 
points pi, P2, pz, • ••, which are poles. 
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The equation (1) is obviously equivalent to the following linear system of 
first order: 

Ui {n -{- r) = — Pi (n) Ui {n -{- r — 1) — Pi (n) Ui (n -\- r — 1) 

— • • • — Pr (n) Ur {n + r — 1) , 



(3) 



U2 {n -{- r) = Ui {n-{- r — 1), 
Us (n -{- r) = U2 {n -{- r — 1), 



Ur {n+ r) = Ur-i {n-\- r — 1 ) . 

In the following paper an analogue of Poincaee's theorem is sought for the 
general linear system: 

(A) yi{n) = aii{n)yi{n — 1) + ■•• + air{n)yr{n— I) (t = 1, 2, ••• , r), 

in which the coefficients a^- (n) have limiting values Aij for n = -{- <^ . This 
extension, which is by no means superficially evident, is contained in Theorems 
1 and 2 of § 2. 

When the coefficients of (A) are constant, it gives for ra = 1 , 2 , • • • the 
indefinite iteration of a collineation. This has led me to the following paper 
on the classification of collineations. 

§1. Preliminary Considerations. 

Consider first the linear substitution (A) for a fixed value of n and seek the 
linear combinations oi yi {n — 1) , • • • , yr {n — 1) which in consequence 
of the substitution are altered only by a multiplicative constant; i. e., 



(4) 



Ylciyi (n) = p^Ciyi (n — 1), 



Inserting here the values of the yt (n) given in (A), we obtain an identity 
which gives the following system of equations for the determination of the 
ratios of the d; 



(5) 



[an (n) — p]ci + Uiiin) C2+ ■ • • + a^ (n) Cr = , 
au (n) ci + [an (n) — p] c^ + ■ • ■ + Ur^ (n) Cr = , 



air (n) ci + a-ir (n) C2 + • • • + [arr (n) — p]cr = 0. 
The elimination of the d gives for p the familiar characteristic equation, 

an(n)—p a^ (n) ■■■ a^i (w) 
an (n) a-n {n) — p ■ ■ ■ art (n) 



(6) 



air (n) a<ir(n) 

Trans. Am. Math. Soc. 28 



fflrr (w) 



= 0. 



344 



E. B. VAN VLECK: ON THE EXTENSION Or 



[July 



When n increases indefinitely, every «„• ( n ) has by hypothesis a Hmit Aij . 
The last equation has therefore the limiting form, 

All — P A21 ••• Arl 

An An — p • • • Ar2 



(7) 



Air 



4 — 



= 0. 



Suppose now that no two roots Ri of (7) are equal, and let them be arranged 
in the order of decreasing modulus, 



(8) 



Rl\^\R2\^ ■■■ ^\Rr 



Then by taking n sufficiently large the r roots of (6) , pi ( ra ) ( i = 1 , • • • , r) , 
may be distinguished from one another by virtue of their proximity to the 
roots Ri of (7). Thus for a sufficiently large value of N we will have 



-Ri — pi{n)\ <\\Ri — Rj 



{n>N;j + i). 



To each p,- (n) there corresponds a set of equations (5), which we will now 
provide with appropriate subscripts and write: 



(9) 



[an (ra) — p,- (ra)] cu (ra) + an (n) ca (ra) + ■ ■ ■ + a^ (n) Cri{n) = 0, 
air (n) cu (n) + a2r (n) Cii (n) + • • ■ + [«„ («) — Pi (n) ]cri (re) = 0. 



From the theory of linear substitutions it is known that r — 1 of these equations 
are independent when the roots of (6) are distinct.* These equations there- 
fore completely determine the ratios of the cy; ( ra ) (j = 1 , ■ ■ • , r) . For a 
like reason the limiting equations. 



(10) 



(^n - Ri) Cu + A21 Cn+ ■■■+ Ari Cri = 0, 

Air Cu +A2rC2i+ ■■■ + { Arr - Ri) Cri = , 



determine the ratios of the Cji . Accordingly for each value of i there is some 
one of the Cjc which may be placed equal to 1 , and the ratios of the other 
Cji {j = 1 , • • • , r) to this particular one will be finite. It follows then, since 
the coefficients of (10) are the limits of those in (9), that these ratios are the 
limits of the corresponding ratios of the cy,- (re) for n = oo . Hence if N is 
chosen sufficiently large, the Cji (n) which correspond to the particular Cji 
taken equal to 1 , may also be set equal to 1. We suppose hereafter that this 
has been done. By means of this convention the Cji (n) (j = 1 , • • • , r) 
will be uniquely determined for re ^ iV, and accordingly we shall have r 

* See, for example, pages 72-73 of Horn's Gewohnliche Dijferentialgleichungen beliebiger 
Ordnung. 
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definite combinations, 

r 

(11) Vi{n) =Y^cn{n)yj{n—\) (i = l,---,r), 

which possess the following property: When transformed by the linear sub- 
stitution (^4), they undergo the multiplicative transformation 

r r 

Hcji{n)yj{n) = pi{n)^Cji{n)yj{n - 1). 
For the limiting form of (11) we have 

(12) Vi{n)=i,CiiYj{n-l) (i = l, .--.r). 

In the theory of linear substitutions* it is shown that the determinant 
I Cji {n) \ can not vanish when the roots of ( 6 ) are distinct. We may therefore 
solve (11) for the y, {n — 1 ) in terms of the Vi{n) . Similarly the determinant 
\Cji\ =1= 0, and equations (12) may be solved for the Yj{n— 1). The 
solutions will be denoted by 

r 

(13) yi{n- \) = Y.dji{n)Vj{n) (i = l, ■••,r), 

(14) r.(n-l) = EZ),,F,(n), 

in which Dy,- is obviously the limit of dji (n) for n = ^ . 

The ratios of the Dji can be most easily determined by considering the linear 
system of difference equations with constant coefficients 

r 

(15) Yi{n) =T.AijYj{n- 1) (i = l, •••,r), 

to which the limiting relations (12) and (14) may be attached. For this 
system we have 

Vi(n+1) = RiVi(n). 
Hence 

(16) Yi(n) = J:DjiVi{n+l) = i: Dji Rj Vj (n) . 

3=1 J=l 

If the expressions for Y,- (ra — 1 ) and Yt (n) which are given in (14) and (16) 
are substituted in (15), and the coefficients of each Vj (n) on the two sides 
of the equation are then equated, we obtain the following r sets of relations 
* Cf. Horn, loc. cit., p. 73. 
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(i= I, ■■■,r): 

{An - Rj) Al + ^12 A-2 + • • • + Air Djr = 0, 
An A-i + (^22 - Rj) Dj,+ ■■■ + A,r T>ir = 0, 

ArX Dyi + .4,2 Z)y2 + • • • + ( ^rr " Rj ) !>,> = . 



(17) 



T^iis shows that the r coefficients Dji (i = 1 , • • • , r) are the homogeneous coor- 
dinates of one of the invariant points of the linear substitution (15). 

In considering the effect of the indefinite increase of n upon any solution 
yi (n) {i = I , ■ ■ • , r) of (A), I propose to deal with them through the 
Vi (n) . For this purpose it will be necessary to find the relations subsisting 
between the Vi (n -\- 1 ) and Vi (n) . With the aid of (A) we obtain 

r r r 

Vi{n+l) = ^Cji{n+l)yj{n) = S Hcji{n+ 1) Ujkin) ykin- 1) 

J=l J=l k=l 

/ ION r r r r 

= S ILcji (re) ajk in)ykin- 1) + J^ Jllcji (re+ 1) 

— Cji {n)]ajk (n) yk (n — 1). 
Now by virtue of (9) we have 

r r r 

Y^^cyi{n)ajk{n) yk{n — 1) = ^ pi{n)cki{n)yk{n— 1) = pi{n)xi{n). 

k=\ J-=1 k=l 

As, furthermore, pi (re) has for re = <» the limit Ri and Cji (re+1 ) — cy» (re.) 
has the limit , equations (18) may be written, 

r 

'Oi{n+ 1 ) = RiVi{n) + ^e'kiyk{n— 1 ) , 

where lim eL ( w ) = for re = <» . Finally, if we substitute for the yu ( n— 1 ) 
their values in terms of the Vj (re) as indicated in (13) and remember that the 
dji have limits for re = co , we obtain, at least for n > N , the system of 
equations, 

r 

(B) Vi {n + 1) = RiVi (n) + '^eij (n) Vjin) (i = l,---,r), 

3=1 

where 

Wvaeij (re) = 0. 

This system of equations forms the basis of the subsequent work. The 
values of the r/,- ( re ) and hence of the Vi{n-\- 1 ) depend upon the initial values 
given to the ^i ( ) and also upon the coefficients aij(p) , ( p = 1 , • • ■ , re ) . 
By exception the yi{n) and Vi{n + 1 ) may all vanish for some values of re , 
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and consequently for all succeeding values. This can happen only if the deter- 
minant of the coefficients an (n) vanishes for at least one value of n . We shall 
assume in § 2 that at least one of the lu {n) does not vanish for an infinity of 
values of n . 

§2. Extension of Poincare's theorem. 
It will be supposed henceforth that 

(19) I ill 1 > i?2 I > • • • > I -Rr 1 . 

When n increases indefinitely, some one of the r moduli | «i (ji) | will be as 
large or larger than the remaining r — 1 for an infinite number of values of n . 
I will show that this \vi{n) | will become and remain dominant over all the 
other r — 1 moduli. 

Suppose first that for an infinite set of values of n we have 

(20) \vi(n) I ^ \vi(n) \ {i = 2, 3, •••, r), 

and denote this set by { n' } . The division of the first of equations (B) by 
»i (n) and of the succeeding equations hy vi {71 -\- 1 ) gives 

A-)1^ ^1 (w+ 1) p , , , 

(21) .a^) =-^' + ^^(^) 
and 

(oo\ Viin+1) _ Rj Vijn) Bj jn) _ _ 

^^-^ vi(n+l) Ri + e^in) vi (n)'^ Ri + e^in) ' > >> 

in which the e ( n ) satisfy for n = n' the inequalities 

(23) \eiin')\ Srein'), \ Bi {n') \ £ re {n') , 

where e ( n' ) denotes the largest of the | etj (n') \m (B). Since lim e ( n' ) = 
for n' = CO , it follows that for a sufficiently large value of n' we have 

\vi{n'+l) \ \Ri\ 
|t,i(n'+l)r = |i?i|+"' 

where a is an arbitrarily prescribed positive quantity subject to the condition 
\Ri\/\Ri\-\-a<. 1 . Consequently n' -\- 1 will be included with n' in the 
set { n' } , provided only n' is taken sufficiently large. It follows that equations 
(20) to (23) will hold without exception from and after some sufficiently large 
value of n. 

From (21) we now obtain 

(24) lim -;—z — = Ri . 
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since lim e (n) = 0._lt may next be shown that 

(25) Hm~-(=0. (i = 2,3, •••,r). 

For if this is not true, there will be at least some one value of i and some 
fixed magnitude M > such that for an infinity of values of n one has 



Call this set of values [ n" } . Then since lim ei (n) = lim et (n) = , 
equation (22) for n = n" may be written 



Viin" + 1) 
v,{n"+l) 






where lim e'i {n") = for n" = oo . Since \ Ri\ <. \Ri\, this shows that 
when n" is sufficiently large, the ratio \vi{n"-\-j) \ / \ Vi (n" -^ j) \ for 
j = 1 , 2 , • • • will decrease in approximately the ratio | R; \ •.{ Ri\ until a 
value of j is reached for which the ratio falls below M . Thus if \vi (n) \ / \Vi (n) | 
is equal to or greater than M for an infinity of values of n , there must also be 
an infinity of such values for which it falls below M . But when it falls below 
M , equation (22) gives 

\viin+ 1)\ ^ \ Ri\ „ , \eiin)\ 



|z>i(n+ 1)1 |fii + ei(n)r^ ' \Ri + ei{n)\ ' 

Since now lim ei{n) = lim e^ (n) = for n = ^ , n may be chosen so large 
in the last inequality that 

\vi{n+l) 



«i(n+l)| 



< M. 



Thus \l \xii{n)\ / I ?)i (n) I < M for a sufficiently large value of n, it must 
remain so for all larger values of n . But this contradicts the hypothesis that 
it exceeds M for an infinity of values of n. Consequently (25) must hold, 
as was to be proved. 

We shall next establish like results in case the r — 1 inequalities 

(26) \v.An)\'^\xi{n)\ (i + 2). 

hold simultaneously for an infinity of values of n . Denote this particular set 
of values by { n' } . Then on dividing equations (2?) for z = 2 by «2 (w) 
and for i =(= 2 by ?)2 ( ^ + 1 ) we obtain the following equations : 

(27) '^^=R. + e.{n), 



1912] A THEOREM OF POINCARfi FOR DIFFERENCE-EQUATIONS 349 

»i(n+l) Ri vi{n.) Bi (n) 



(28) 



1)2 ( n + 1 ) R2 + 62 (n) 1)2 ( n ) -R2 + 62 ( n ) ' 



(29) — -, rTT="5 — t 7~T ^^+75 — 1 T^ (t = 3, 4, • • • , r), 

V2in+l) -R2+ 62(71) »2(n) ii2+ 62(71) 

in which the e (n) f or n = n' satisfy the inequaHties 

(30) \ei{n')\Sre{n') (i = 1, 2, ••• , r), 

where e {n') denotes the largest of the | 6,7 («,')[ in (B). 

Equations (29) show that for a sufBciently large value of n' we have 

(31) \vi{n'+l)\<\v2{n'+l)\ (i=3,4, ....r) 
since 

\Ri\<\R2\, \viin')\ £\v2{n')\, lime(7i') = 0. 

n^—v) 

Furthermore, we must also have for a sufBciently large value of n' 

(32) \vx(n'+l)\<\v2(n'+l)\, 
for otherwise the inequality 

hi(n'+l)| ^\vi{n'+l) I (i = l,3, •••,r) 

would hold for an infinity of values of n' . But the latter is impossible, since 
it has been shown that \vi{n) \ would then ultimately become and remain 
dominant over all the other \vi{n) | , which contradicts the hypothesis that 
(26) is to hold for an infinity of values of n . Now (31) and (32) together show 
that n' + 1 will be included simultaneously with n' in the set { n' ] , provided 
only n' is sufficiently large. It follows that (26)-(30) will all hold without 
exception from and after some fixed but large enough value of n . 
From (27) we see next that 

(33) lim ~, — r — = R2 . 

From (29) it follows also that 

(34) lim ''i''t !! =0 (i=3, ••-.r), 

n=» ■»2 (m+ 1) 

the reasoning by which this is established being exactly parallel to that by 
which (25) was established on our former hypothesis. Lastly, it remains 
to be proved that 

(35) lim^=0. 

Suppose, if possible, this not to be true. Then there must be some magnitude 
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M > which \ vi (n) \/ \v2 (n) \ will exceed for an infinity of values of n . 
Denote this set of values by { n" } . For these values (28) may be written 



v,in"+l) \R^^' ^"^ V 



vAn")' 



where lim e' (n") = for n" = ^ . Since | i^i | > | i^2 | , this shows that 
ultimately n" + 1 will be included with n" in the set { n" ] . Then as i 
increases through the natural number system, the ratio |i'i(ji"+t)| / \vi{n"+i\) 
will increase successively by factors approximating closer and closer to 
I i?i I / I i?2 I , and hence finally we will have \ Vi (n" + i) \ > \vi {n" + i) \ 
from and after some fixed value of i. But this contradicts our hypothesis 
(26). Hence (35) must hold, as stated. 

In similar manner we may treat successively f or j = 3 , 4 , • • ■ , r the hy- 
pothesis that the r — 1 inequalities 

\vj{n)\^\vi{n)\ (i+i) 

hold for an infinity of values of n . On dividing (B) by Vj (n) and Vj (n-'r 1) 
we obtain three sets of inequalities similar to (27), (28), (29) , which furnish 
like results. We reach thus the following conclusion: 

Theorem 1. If in the system of differential equations (B) we have 

lim eij (n) = {i, j = 1, 2, ■■■ , r) 

W=co 

and we suppose 

\Ri\>\R,\> ■■■ >\Rr\, 

then as n increases indefinitely, either by exception the r quantities Vi{n) loill all 
vanish from and after some fixed value of n, or some one of them will ultimately 
dominate over the other r — 1 . If this he Vk (n) , then 

Vk+iin) 
lim -. — ^ = Rk 

»=« Vk{n) 

liin^_d_^) (i + fc). 

M==o Vk{n) 



and 



We return now to our initial equations (A). We have seen that in the re- 
lations 

r 

(13) yiin-1) = 'E,dji{n)vj{n) (i = l,---,r) 

the coefficients tend toward limiting values Dji . The determinant ] Dji | 
can not vanish since it is the reciprocal of ] Cji ] 4= . Consequently for 
each value of j not all of the Dji {i = 1 , • • • , r) can vanish. Consider now 



1912] A THEOREM OF POINCAE6 FOR DIFFERENCE-EQUATIONS 351 

the equation 



1^dji(n+ l)vj{n+ 1) 
(36) 



Viin) j=i 



yi{n- 1) ^ 



If, as in Theorem 1, Vk (n) denotes that one of the Vi (n) which ultimately 
becomes dominant, and if Dki denotes one of the Dkj (j = 1 , • ■ ■ , r) which 
does not vanish, it follows by application of the theorem to (36) that 

yi{n—l) dki{n+l)vk{n+l) »/, (n+1) 

lim — ,— = hm T — -. — r — 1-^ = lim -r—. — = lik ■ 

»=« yi{n) dki{n)Vk{n) Vk{n) 

Also we have 

yi{n) dki{n+ l)vk{n+ \) Dki /•,,^ 

lim — ;^--=hmi — -, j— TT -, j— rT=n-- (»4=0. 

This shows that the point yt (n) (i = 1 , • • • , r ) for w = co approaches as 
its limit one of the invariant points of (15). This result may be restated as 
follows : 

Theorem 2. Let (A) be a systevi of difference equations in which the coef- 
ficients an (n) have limiting values A a for n= <x> . If the roots Ri of the limit- 
ing characteristic equation (7 ) are no two of equal modulus, then unless by exception 
the yt (n) all vanish from and after some fixed value of n , there is at least one 
function yi (n) in the solution yj (n) , (j = 1 , • • • , r) , for which we have 

y yi{n) 

hm — ;; r^ = Rk , 

„=^yi (w — 1) 

where Rk is one of the roots of (7). 

If also {A) is a sequence of non-singular collineations in r homogeneous vari- 
ables, tvhich have for n = <^ a limiting form, then any initial point y\ ( ) , 
• • • yr (0) by the successive application of the collineations will approach as a 
limit one of the invariant points of the limiting collineation (15). 

The theorem of Poincare is an immediate consequence of the first part of 
this theorem, for by the very construction of (3) we have 

,. Uj (ra) Ui-ijn) 

hm — -. — r-^r\ = hm -. z-r {^ = 2, ■■• ,r). 

«=« Ui{n+ 1 ) M,-_i ( n — 1 ) 

To guard against misconception I will add in conclusion that not all of the 
r ratios yt (n) / y^ (n — 1 ) need have limits f or n = co , and even if they 
have, not all of the limits need to be roots of (7). The truth of the latter 
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statement will be obvious from a consideration of the special system, 
Di (n) = RiVi{n— 1), 
V2 {n) = e {n) vi (n — I) + R2V2 {n — 1) 

in which e (n) is an infinitesimal having the value 

(Si- E2) !'i^~ ]] (0<R,<S,<Ri). 

^\{n — 1 ) 

Here the limit oi V2 (n) f V2 {n — 1 ) for the particular integral ».• ( n ) = » j (n) 
is an arbitrarily chosen constant »Si . 
Univehsitt of Wisconsin. 



